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Abstract. I cross-correlate the OGLE database of 
Cepheids in the Large and Small Magellanic Clouds (MCs) 
with the second incremental release of the 2MASS JHK s 
survey, and the denis IJK s point source catalog towards 
the MCs. 

Of the 3384 Cepheids in OGLE in both Clouds, 1745 
have a counterpart in the 2MASS survey within an 1" ra- 
dius and good-photometry in all three bands. Only 173 
have a counterpart in the DENIS survey within an 1" ra- 
dius and good-photometry in all three bands. The reason 
for this difference is that the limiting magnitudes of the 
2MASS survey are considerably fainter than for the DENIS 
survey. 

The standard stars of the Carter IR photometric sys- 
tem are also correlated with the 2MASS survey to derive 
transformation equations from the natural 2MASS system 
to the Carter-system. 

In order to describe to first order the geometrical con- 
figuration of the MCs a plane is fitted to the reddening- 
free Wesenheit index (W = I - 1.55 x (V — I)) and the 
inclination and position angle of the line-of-nodes are de- 
termined. For the LMC an inclination angle of 18 ± 3° is 
derived, somewhat smaller that previous determinations. 
For the SMC a value of i — 68 ± 2° is derived in good 
agreement with previous determinations. These results are 
being used to take out the effect that some Cepheids are 
closer to the observer than others. For the LMC this ef- 
fect is small over the area covered by the OGLE survey, 
but for the SMC PL-relations with a reduced scatter are 
derived in this way. Nevertheless, the dispersion in the 
PL-relation in W for the SMC is larger than that for the 
LMC, indicating a significant intrinsic depth of the SMC 
of about 14 kpc (front-to-back). The depth of the LMC is 
found to be small compared to the intrinsic scatter in the 
PL-relation. 

Using single-epoch 2MASS data for the OGLE Cepheids 
PL-relations are fitted in JHK (on the Carter system) for 
both Fundamental mode (FU) and First-overtone (FO) 
pulsators, taking out the effect of the inclined disk. Be- 



cause of the large number of stars available, this results in 
the most accurately determined slopes in the infrared up 
to now, and slopes for FO pulsators derived for the first 
time. This provides additional constraints for theoretical 
models. 

The 16 derived PL-relations (in WJHK for FU and 
FO pulsators for LMC and SMC) are used to derive a 
relative distance modulus Asmc-lmc = 0.50 ± 0.02 as- 
suming no metallicity correction. Two recent theoretical 
models predict different metallicity corrections, indicating 
a systematic uncertainty of up to 0.1 mag. 

Using the database of Galactic Cepheids observed 
with hipparcos and the method outlined in Groenewe- 
gen & Oudmaijer (2000), the zero points of Galactic PL- 
relations for FU pulsators in W and K are derived, using 
the newly derived slopes for LMC and SMC pulsators. 

Combining the zero points of the Galactic and MCs 
PL-relations distance moduli of 18.60 ± 0.11 (based on 
W) and 18.55 ± 0.17 (based on K) to the LMC, and of 
19.11 ± 0.11 (W) and 19.04 ± 0.17 (K) to the SMC are 
derived (without taking into account possible metallicity 
corrections) . 
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1. Introduction 

Cepheids are important standard candles in determin- 
ing the extra-galactic distance scale. The results of the 
HIPPARCOS mission allowed a calibration of the Galactic 
Period-Luminosity (PL—) relation in the U-band (Feast & 
Catchpole 1997), V and /-bands (Lanoix et al. 1999a) and 
V, I, K and the reddening-frec Wesenheit index (Groe- 
newegen & Oudmaijer 2000, hereafter GO00). Feast & 
Catchpole (1997) and GO00 combined these data with 
available data on LMC Cepheids to determine the dis- 
tance to the LMC. It is well know that the infrared has 
some considerable advantages over the optical region for 
deriving distances. First of all the problem of reddening is 
less important, and second, as summarised in GO00, the 
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effect of metallicity on the PL-relation seems to be less. 
Of the 236 Cepheids in the sample considered in GO00, 
only 63 have sufficiently well determined infrared light 
curves to determine accurate intensity-mean magnitudes 
(also see Groenewegen f999). In addition, the number of 
Magellanic Clouds (MCs) Cepheids with well determined 
infrared lightcurves is also rather small, of order 20 in each 
Cloud (Laney & Stobie, 1994, 1986a). 

The micro lensing surveys have increased the number 
of known Cepheids in the MCs tremendously and they 
provide optical photometry for many of them. In particu- 
larly, Udalski et al. (1999b, c; hereafter U99b,cQ) recently 
published the results on Cepheids in the Small and Large 
Magellanic Cloud from the OGLE survey. Complementary, 
in the infrared, the 2MASS survey is an ongoing all-sky 
survey in the JHK S near- infrared bands, that, when com- 
pleted, will contain of order 300 million stars (Beichman 
et al. 1998). The 2MASS team released their second incre- 
mental data set in March 2000, that includes the MC area. 
In parallel, the DENIS survey is a survey of the southern 
hemisphere in IJK S (Epchtein et al. 1999), that recently 
released a Point Source Catalog (PSC) of sources in the 
direction of the MCs (Cioni et al. 2000a). 

In this paper the OGLE database is correlated with the 
2MASS second incremental database and the denis PSC of 
MC stars (Sect. 2), transformation formulae from the nat- 
ural 2MASS system to the Carter (1990) system are derived 
(Sect. 3), the spatial structure of the MCs is investigated 
(Sect. 4), and PL-relations in JHK for fundamental mode 
and first overtone pulsators derived (Sect. 5). The results 
are discussed in Sect. 6. 



2. OGLE, 2MASS and denis Sample selection 

U99b,c describe the dataset regarding Cepheids in the di- 
rection of the LMC and SMC, obtained in course of the 
OGLE-n micro lensing survey. Twenty-one fields in the 
central parts of the LMC, and 11 fields in the central 
parts of the SMC of size 14.2' x 57' each were observed 
in BVI, with an absolute photometric accuracy of 0.01- 
0.02 mag. U99b,c present tables per field, with (amongst 
other items) identification number, pulsation period, BVI 
photometry, E(B—V) and pulsation mode. The E(B — V) 
value was derived by them using red clump stars, and the 
pulsation mode was derived by them from two Fourier pa- 
rameters derived from the /-band lightcurve (see U99b,c 
for details). These data were obtained from the OGLE 
homepage ( http://sirius.astrouw.cdu.pl/~oglc/ ) in April 
2000. 

As noted in U99b,c there are about a 100 objects in 
each cloud that appear in 2 or more fields. They provide 



1 In April 2000 the OGLE team released a re-reduced dataset 
changing the photometry slightly, and hence the PL-relations 
with respect to those in the papers published in 1999. The new 
data and corresponding PL-relations can be found on the OGLE 
Homepage. When appropriate I will refer to these new values. 



a list with the relevant field and identification number, 
but I did this positional correlation independently. Values 
for the photometry, period and E(B — V) were averaged 
where appropriate. The final number of Cepheids detected 
by OGLE in the LMC is 1335, and in the SMC 2048. 

The 2MASS survey is an ongoing single-epoch all-sky 
survey in the JHK S near-infrared bands. On March 2, 
2000 the 2MASS team released the second incremental data 
set, that includes the MC area. The easiest way to check 
if a star is included in this release is by uplinking a source 
table with coordinates to the 2MASS homepage. Such a ta- 
ble was prepared including the 1335 LMC and 2048 SMC 
Cepheids in OGLE and 236 Cepheids in the HIPPARCOS 
catalog considered in GO00. 

The correlation between the source list and the 2MASS 
database was at first instance done using a search radius 
of 5". Inspection of the match between the OGLE and HIP- 
PARCOS Cepheids and 2MASS indicated a much better po- 
sitional coincidence. This is to be expected as the 2MASS 
positional accuracy is <0.2" (Cutri et al. 2000) and the in- 
ternal accuracy of the OGLE database is about 0.15" with 
possible systematic errors of <0.7 " (U99b,c). In quite a 
few cases there was a second 2MASS source (not the coun- 
terpart) within the search radius of 5". The match between 
2MASS and HIPPARCOS was excellent with the largest posi- 
tional difference being 0.56". In view of this, a second and 
final correlation was made for the HIPPARCOS and OGLE 
Cepheids using a search radius of 1". Seventy- two correla- 
tions were found with the HIPPARCOS Cepheids, 1511 with 
the OGLE SMC, and 894 with the OGLE LMC Cepheids. 
The small number of correlations with HIPPARCOS is due 
to the fact that the 2MASS survey is not yet complete in 
sky coverage. 

On this 2MASS dataset a further selection was per- 
formed to retain only those sources with reliable JHK 
photometry in all three bands. It is worthwhile to recall 
here that the "rd-flg" in the 2MASS database gives essen- 
tial information on the quality of the magnitudes (Cutri 
et al. 2000). In particular, objects where the rd-flg is 
(indicating an upper limit) or 3 or 8 (indicating satura- 
tion) in any of the three bands are eliminated. The fi- 
nal 2mass database of Cepheids consists of 52 hippar- 
COS Cepheids (not discussed further), 825 OGLE LMC and 
920 OGLE SMC Cepheids. The 2MASS and OGLE coverage 
of the LMC are not identical, and no 2MASS sources are 
found in OGLE fields LMC-SC6, 7 and 21. 

The denis survey is a survey of the southern sky 
in IJK S (Epchtein et al. 1999). In April 2000, a point 
source catalog was released containing sources and single- 
epoch photometry in the direction of the MCs (Cioni et al. 
2000a). A subset of 264347 LMC and 69829 SMC sources 
detected in all three band was obtained using FTP from 
www.strw.leidenuniv.nl (cd /pub/ldac/dcmc). These were 
correlated against the OGLE survey on position. As both 
OGLE and DENIS contain /-band data it was possible to 
verify how the average difference between the OGLE and 
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denis /-data depends on the search radius. Since the OGLE 
data contain intensity-mean magnitudes, but the denis 
data are single-epoch data, one has to allow for variabil- 
ity. For the SMC, and a search radius of 2", 93 cross- 
correlations were found. The average difference in posi- 
tion is 0.68" and the mean absolute difference in I is 0.18 
mag. However, for the 21 stars with a positional difference 
larger than 1", 9 (or 43%) have a difference in / of larger 
than 0.8 mag, which is a typical value for the amplitude for 
a Cepheid in the /-band. For the stars with a positional 
differences less than 1", this is only 5 out of 67 or 7%. 
Given the internal accuracy of the positions in the denis 
survey of 0.5" (Epchtein et al. 1999), a search radius of 1" 
seems appropriate to reduce the number of spurious asso- 
ciations, which is supported by the large fraction of highly 
deviant /-band data when a larger search radius is used. 
As discussed below some true identifications are missed 
in this way however. Also some false coincidences are un- 
doubtedly included (as for the correlation of OGLE with 
2mass) but these will be removed using sigma-clipping 
when deriving PL— relations. Using a search radius of 1" 
the number of correlations between denis and OGLE is 71 
in the SMC, and 102 in the LMC, for stars with good 
photometry in all three bands 

It was verified whether all correlations between denis 
and OGLE were also found in the match between 2MASS 
and OGLE. Of the 71 SMC OGLE + denis stars, 61 are also 
found in 2MASS. In three of the 10 cases where there is no 
match this was due to the positional difference which was 
> 1" between OGLE and 2MASS. In the other 7 cases it 
is less clear, possibly related to the variable nature of the 
objects. Of the 102 LMC ogle + denis stars, 80 were 
also found in 2MASS. In one of the 22 other cases this 
was due to a positional difference of > 1" between OGLE 
and 2MASS, and < 1" between OGLE and DENIS. One case 
is puzzling, and the other 20 cases are all in OGLE fields 
LMC-SC6, 7, 21, which apparently have not been covered 
by 2 MASS. 

The number of correlations beween DENIS and OGLE is 
much less than the number of correlations between OGLE 
and 2MASS. Some tests using denis data of sources de- 
tected in / and J only, showed that this is mainly due to 
the lack of detections in K. Some tests in a few OGLE fields 
showed that there may be an equal number of Cepheids 
in OGLE with DENIS / and J photometry. In any case it 
is clear however, that for the purpose of determining PL- 
relations in the infrared using the largest possible samples, 
the 2MASS survey is of more interest, because it reaches 
fainter magnitudes. 

In Appendix A a comparison is presented between the 
DENIS / and OGLE / and denis JK and 2MASS JK pho- 
tometry. This has its limitations since the comparison is 
done using variable stars. Yet it may still be of some inter- 
est to other workers, as such a comparison has not been 
done yet. 



3. Transformation of the photometric system 

The photometry presented by the 2MASS-team is on the 
natural system. However, many of the extisting infrared 
photometry on both Galactic (Laney & Stobie 1992, Groe- 
newegen 1999) and MC Cepheids (Laney & Stobie 1994) 
is on the Carter (1990) system or has been transformed 
to that system. It is therefere imperative to estimate any 
systematic differences between the Carter and the natural 
2MASS system. 

From Carter (1990) and Carter & Meadows (1995) 116 
stars with K > 5 were taken, and a table with coordi- 
nates was uplinked to the 2MASS homepage to perform 
a search on position using a search radius of 2". Thirty- 
four matches were found; the largest positional difference 
was 1.2". However, not all entries are useful as some of 
the brighter stars still saturate the 2MASS detectors in 
one or more bands, despite the a-priori exclusion of stars 
with K < 5. The vast majority of stars in 2MASS have 
"default magnitudes obtained using profile-fit photometry 
performed simultaneously on the combination of all six in- 
dividual 1.3-seconds "Read 2-Read 1" (R2-R1) exposures" 
(Cutri et al. 2000). These sources have rd-flg = 2 in the 
appropriate band. Sources brighter than 7-8 magnitudes 
will saturate in the 1.3-s R2-R1 exposures. These objects 
have default magnitudes from aperture photometry per- 
formed on the 51-millisec Rl frames. Such sources have 
rd-flg = 1 in the appropriate band. Stars brighter than 
approximately fifth magnitude will saturate in even the 
51-millisec exposures. The rd-flg value for the appropriate 
band is 3. 

Among the 34 sources, there are 15 sources with rd-flg 
= 111, 10 with rd-flg = 112, 1 with rd-flg = 122, 4 with 
rd-flg = 222 and 4 sources with rd-flg = 3 in atlcast one 
of the bands. The latter were not considered. 

Figure 1 plots the differences "Carter - 2Mass" versus 
(J - K) (2Mass), and "2Mass - Carter" versus (J - K) 
(Carter) in JHK. For these nearby stars reddening is ne- 
glected. The errors in the individual magnitudes are taken 
from the 2MASS database (typically between 0.01 and 0.03 
mag), or arc 0.025 mag in the case of the Carter system fol- 
lowing the estimate in Carter & Meadows (1995). Linear 
least-square fits show that the significance of the derived 
slopes is less than la in HK and about 1.5er in K, indi- 
cating that there is at present no indication of a colour 
term in the transformation from the 2MASS to the Carter 
system. The average off-sets are (Carter) — K (2Mass) 
= +0.017 ± 0.008, //(Carter) - //(2Mass) = -0.021 ± 
0.009 and J(Carter) - J(2Mass) = +0.048 ± 0.007 (in- 
ternal errors). The dispersion in all three plots is about 
0.03 mag. This analysis indicates (A) that there might be 
small off-sets between the two infrared systems, and (B) 
that the absolute photometric calibration of both is uncer- 
tain at the 0.01-0.02 mag level. This issue may addressed 
in more detail when the 2MASS survey has been completed 
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Fig. 1. Comparison between JHK in the Carter (1990) and the natural system of 2MASS. 



and this analysis may be repeated using a larger number 
of stars. 

Regarding the same issue, Nikolaev & Weinberg (2000) 
note that the 2MASS system is similar to the CIT/CTIO 
system (Elias et al. 1982, Persson et al. 1998) except for 
the K s band. They note that the difference of the absolute 
value of (ifciT — K s ) is less than 0.05 mag. 

4. The spatial structure of the LMC and SMC 

In this section the spatial structure of the MCs is being 
investigated in terms of the usual approximation of a flat 
disk, using the largest possible dataset, that is the OGLE 
sample (without the correlation with 2mass). The sec- 
ondary aim is to correct for any depth effects in order to 
obtain tighter PL-relations. 

Cepheids have previously been used to determine the 
spatial structure of the MCs, for example Caldwell & 
Coulson (1986) whose model for the LMC is a plane of 



insignificant thickness inclined by 29 ± 6° with the closest 
part at position angle (P. A.) 52 ± 8°. For the SMC they 
find that a planar model is barely sufficient to describe this 
galaxy, and derive an inclination of 70 ± 3° with the clos- 
est part at P. A. 58 ± 10°. Laney & Stobie (1986a) reach 
similar conclusions and derive i = 45 ± 7 for the LMC 
with the nearest part at P. A. 55 ± 17°. Values for inclina- 
tion and position angle derived from other methods and 
a more general discussion can be found in the monograph 
by Westerlund (1997). Most recently, Weinberg & Niko- 
laev (2000) used star counts from the 2MASS survey of the 
LMC, and a model for the number density of stars, to de- 
rive inclinations between 24 and 28 degrees, and position 
angles of the line-of-nodes between 169 and 173 degrees 
depending on the population, and an inclination of 42.3 
± 7.2° from carbon-rich Long Period Variables (LPVs). 

The coordinate transformations that are necessary 
have conveniently been written down in Appendix A.l 
in Weinberg & Nikolaev (2000). To summarise, for every 
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Cepheid under consideration one knows its right ascension 
and declination, and one can derive its distance r (see be- 
low). This coordinate system (a, S, r) is transformed into 
a rectangular coordinate system (xo, yo, zq), which has its 
origin at the center (ao, So, R) of the LMC or SMC, the 
zo-axis towards the observer, the xo-axis anti-parallel to 
the right ascension axis, and yo parallel to the declination 
axis. Another rectangular coordinate system is introduced 
(x 1 , y', z'), which is rotated about the zo-axis by an angle 
9 counterclockwise (counting East of North), and about 
the new x'-eods (the line-of-nodes) by the inclination an- 
gle i clockwise (that is rotation from the new y-axis to- 
wards the —zo-axis). Note that there seem to exist slightly 
different definitions in the literature regarding this issue, 
related to, for example, the definition of negative/positive 
inclination^} 

The above procedure assumes one knows the distance 
to the center of the galaxy, and the distance to the individ- 
ual Cepheids. To this end an observed PL-relation is used, 
together with a typical distance (the assumed value for R 
is arbitrary because one basically look at residuals in mag- 
nitude space relative to the observed mean PL-relation. 
However, for convenience, it is preferred here to work in 
physical, rather than magnitude, space, and therefore an 
arbitrary mean distance to the center of the galaxy has to 
be adopted). 

In particular, the PL-relation using the reddening-free 
Wesenheit-index using the V and / colours is being used 
as the reference PL-relation. This is done as the observed 
scatter in this relation is less than in the PL-relations in 
V or /, and because the Wesenheit-index, by construction, 
does not depend on the assumed value for E(B — V). It 
depends only on the reddening law. Its definition is W 
= I - 1.55 x (V — I) for the OGLE filter system and the 
extinction law by Cardelli et al. (1989) (Udalski et al. 
1999a; hereafter U99a) 

4.1. The LMC 

U99a derived the zero point and slope of the PL-relation 
in W, for periods longer than about 2.5 days and funda- 
mental pulsators (FU), and applying clipping at the 2.5a 
level (see Table 1, entry 1). They used this cut-off in period 
as in the SMC there is a change of slope at shorter peri- 
ods (Bauer et al. 1999). Using the same criteria I derive 
almost identical results (solution 2), altough I would pre- 
fer to use a less stringent clipping to have a larger sample. 
A solution with clipping at the 4cr level is also included 
(solution 3). 

The change of slope observed in the EROS data for the 
SMC was not observed in the LMC (Bauer et al. 1999), 

2 In particular, using the present notation, the result of Cald- 
well & Coulson (1986) implies 8 = 232 ±8, the result of Laney 
& Stobie (1986a) implies 6 = 235 ± 17, and that of Weinberg 
& Nikolaev (2000) implies a P.A. of the line-of-nodes between 
259 and 263 degrees for the LMC. 



and also theoretically, one would expect this change of 
slope to occur at shorter periods for higher metallicity (Al- 
ibert et al. 1999). Table 1 also includes a fit using no con- 
straint on the period (solution 4). The results are nearly 
identical, but obtained for a larger sample. 

U99a only briefly mention the result for the slope de- 
rived for First Overtone (FO) pulsators, and, in fact, more 
details are given on the OGLE Homepage (entry 5 in Ta- 
ble 1), which has been derived for stars with logP > —0.2 
(Udalsky, private communication). I find essentially the 
same result (solution 6), also using a less stringent 4cr 
clipping (solution 7). However, for consistency, a cut-off 
at logP = 0.4 for FU pulsators should correspond to a 
cut-off at about logP = 0.25 for FO pulsators, according 
the usual ratio of Po/Pi (for example Feast & Catchpole 
1997). Table 1 gives the solutions for stars above and be- 
low this cut-off in period (solutions 8, 9). The slopes are 
nearly identical, but there is a significant off-set in zero 
point. Bauer et al. (1999) for a smaller sample of 113 FO 
LMC Cepheids detected by EROS found that both zero 
point and slope did not depend on a cut-off in period. 
However their error bars are roughly 4 times larger than 
here, and so their sample may not have been large enough 
to detect this effect. 

For the investigation of the effect of an inclined disk 
the FU pulsators are being considered first, using the 
relevant observed PL- relation (solution 4). An arbitrary 
distance to the center of the LMC of R = 50 kpc is 
assumed, and the center of the disk is assumed to be 
given by the mean of the 713 Cepheids in the sample: 
a = 81.384, S = -69.782 degrees (J2000). The observed 
PL-relation is then transformed into an absolute relation 
M\v = —3.320 logP —2.615, and with this the distance to 
every Cepheid can be calculated, and the coordinates (xq, 
yo, Zq) and {x 1 , y' , z') determined for an assumed position 
angle and inclination. 

There may be one caveat of working in real space 
rather than magnitude space that should be pointed out 
here. The dispersion around the observed PL-relation is 
due to three effects. The depth effect along the line-of- 
sight, or, equivalently, the distance zq to the plane of the 
sky, photometric errors, and "cosmic" scatter. The latter 
term arises because a linear PL-relation is an approxima- 
tion, metallicity effects, or other effects. The photometric 
errors are small (of order 0.01 mag) and the other two 
effects are difficult to distinguish, and when considering 
correlations of the residual magnitudes with respect to 
the mean PL-relation with position need not to be dis- 
tinguished. When working in physical space, the z val- 
ues derived are NOT the actual distances below or above 
the mid-plane. In fact the values for the distance to the 
mid-plane are upper limits as the "cosmic" dispersion is 
non-zero. 

Figure 2 shows the distributions for i = 6 = 0. The 
upper- left panel shows in filled squares stars that are closer 
to the observer (that is brighter than the PL-relation) and 
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Table 1. 


PL-relations in the LMC of the form 


M = a 


x logP 


+ b 


solution 


M 


a 


b 


a 


N 


Remarks 


1 


W 


-3.300 ±0.011 


15.868 ± 0.008 


0.058 


668 


U99a; FU; log P > 0.4; 2.5cr clipping 


2 


W 


-3.303 ±0.012 


15.869 ± 0.008 


0.059 


639 


This paper; FU; log P > 0.4; 2.5cr clipping 


3 


W 


-3.310 ±0.014 


15.871 ± 0.009 


0.073 


673 


This paper; FU; log P > 0.4; 4a clipping 


4 


w 


-3.320 ±0.013 


15.880 ± 0.009 


0.074 


713 


This paper; FU; all log P; 4<j clipping 


5 


w 


-3.425 ±0.017 


15.380 ± 0.006 


0.060 


457 


U99a; FO; logP > -0.2; 2.5cr clipping 


6 


w 


-3.434 ± 0.017 


15.383 ± 0.006 


0.060 


424 


This paper; FO; logP > —0.2; 2.5cr clipping 


7 


w 


-3.400 ± 0.018 


15.369 ± 0.007 


0.072 


450 


This paper; FO; all logP; 4a clipping 


8 


w 


-3.300 ± 0.032 


15.323 ± 0.013 


0.068 


344 


This paper; FO; log P > 0.25; 4cr clipping 


9 


w 


-3.266 ± 0.050 


15.385 ± 0.009 


0.084 


108 


This paper; FO; log P < 0.25; 4a clipping 


10 


w 


-3.337 ±0.013 


15.890 ± 0.008 


0.072 


713 


FU; all log P; 4a clipping; planar effect taken out 


11 


w 


-3.398 ±0.018 


15.370 ± 0.007 


0.072 


450 


FO; all logP; 4a clipping; planar effect taken out 



in open circles Cepheids that are further away (fainter 
than the PL-relation)^| The other 3 panels show the pro- 
jection on to different planes. From the top right panel a 
rotation angle of approximately 26 (c.q. 206) or 116 (c.q. 
296) degrees is derived. The question of the inclination 
is not immediately clear. Figure 3 shows the distribution 
and projections for 9 = 26 and i = 0. The x'-axis now 
runs along the major axis from bottom left to top right 
and y'-axis runs along the minor axis from bottom right to 
top left. The x' and y' coordinates are now uncorrelated, 
by construction. Their is no correlation between z' and 
y' , indicating no rotation over the x'-axis. The correlation 
between z' and x' indicates a rotation along the y'-axis 
by 18 degrees with the south-east part tilted towards us. 
Given the definition of the coordinate system, and consid- 
ering the errors in the slopes, this implies a position angle 
of the line-of-nodes of 9 = 296 ± 0.5, i = 18 ± 3°. 

This analysis was repeated using all FO pulsators (so- 
lution 7), and Fig. 4 shows the result for i — and 9 = 24. 
The inclination implied by the slope in the z' and y' plane 
is not significant being 16 ± 17°, and the slope in the other 
plane of 10 ± 5° is consistent with that derived from the 
FU pulsators. 

These results are in fair agreement with previous re- 
sults in the literature, where larger inclinations are usually 
quoted (see the beginning of Sect. 4.0). Probably the main 
reason is that those studies covered a much larger area, 
and not only the bar region. In particular it would be ex- 
tremely useful to have OGLE fields at some distance from 
the bar, along P.A. of about 30, 120, 210 and 300 degrees. 
One other possibility is that the structure of the LMC is 
truly different in the central region from the more global 
structure. A very recent preprint (Zhao & Evans 2000) 
eluded to the possibility of an off-setted and misaligned 
LMC bar with respect to the LMC disk. 

As a further illustration Fig. 5 depicts the distance 
from the mid-plane assuming no intrinsic scatter of the 
PL-relation. In fact, this plot illustrates the upper limit 



3 In the electronic version, this is highlighted by using differ- 
ent colour schemes as well. 




-10 1 -10 1 

Y' (kpc) X' (kpc) 



Fig. 2. Distribution on the sky of FU Cepheids who are 
brighter (filled squares) or fainter (open circles) than the 
mean PL- relation (upper left panel), and projection in 
(x',y',z') coordinates for i = and 9 = 0. Solid lines 
indicate the results of linear least-square fits. 



to this distribution, and the true distribution (x-scale) 
should be compressed by a factor of about 9 (see below) . 

From McCall (1993) it is estimated that the layer con- 
taining 99% of young stars has a thickness of 1.0 ± 0.3 
kpc, for a distance to SN 1987A of 51.2 ± 1.2 kpc (Pana- 
gia 2000). The half-thickness of the layer containing 50% 
of the young stars can then be calculated to be 0.1 ± 0.03 
kpc. McCall also derives H/R = 0.120 ± 0.031, where H 
is the vertical scale height and Rq the radial scale length 
of the LMC disk. For Rq between 1.4 and 1.8 kpc derived 
from 2MASS data (Weinberg & Nikolaev 2000), this results 
in H = 0.19 ±0.06 kpc. Taking this latter value, it is con- 
cluded that at the distance of the LMC the thickness of 
the LMC disk gives rise to a la dispersion in PL-relations 
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Fig. 3. As Fig. 2 for i = and = 26°. 
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Fig. 4. As Fig. 2 for i = and 6* = 24° and FO pulsators. 

of approximately 0.008 mag only. The total front-to-back 
depth corresponds only to 0.043 mag. 

A final remark regarding this topic is that the the- 
oretically predicted dispersion in the Wesenhcit-index is 
0.10 mag (Bono & Marconi, 2000, private communica- 
tion). This is in fact slightly larger than the observed one. 
However, this is again consistent with the notion that the 
intrinsic depth of the LMC plays no role in the observed 
scatter in the PL- relations towards the LMC. 

Weinberg & Nikolaev (2000) recently claimed the 
LMC to be geometrically thick along the line of sight with 



-4 -2 2 4 

Z (kpc) 

Fig. 5. If the intrinsic spread of the PL-relation were zero, 
this would be the distribution of stars along the line-of- 
sight towards the LMC. Positive- z are towards the ob- 
server. As the intrinsic spread is in fact 0.008 mag, the 
horizontal scale should be compressed in such a way that 
the 4 kpc thick actually correponds to 450 pc. 



a thickness of ~14 kpc. They used stars in a color selected 
region based on 2MASS data and assumed these to be 
carbon-rich LP Vs. They derive widths of the magnitude 
distribution functions of 2a = 0.7 mag, while they claim 
their standard candles to have an intrinsic spread of 0.2 
mag, from which they derive a front-to-back thickness of 
~14 kpc, after they rejected other possible explanations. 
For a similar population, Bessell et al. (1986), from the 
intrinsic line-of-sight velocity dispersion of old LPVs, 
derived a vertical scale height H ~0.3 kpc, consistent 
with the notion of a thin disk. 

Although it will turn out academic for the LMC, it 
will be discussed now, how an improved PL-relation can 
be derived once position angle and inclination have been 
determined. If the approximation of a thin disk is correct, 
one can calculate those distances, r', respectively, r, that 
make zq = and z' = (for fixed a and 8 and given 9 and 
i). That is, for every star one can subtract the correction 
in magnitude, 5 log(f/r'), to put all stars in the plane of 
the sky, and then fit the PL-relation again. Planar effects 
are in this way taken out. The result is listed in Table 1 
as solutions 9 and 10, and are essentially the same as so- 
lutions 4 and 7 with only marginally lower errors and the 
fit coefficients. This indicates that over the region covered 
by OGLE, depth effects are not important, and that the 
width of the PL- relation is intrinsic. The resulting final 
PL-relations are shown in Figs. 6 and 7. 

4.2. SMC 

This section is largely a repeat of the previous subsection 
but for the SMC. U99a derived the zero point and slope of 
the PL-relation in W, for periods longer than about 2.5 
days and fundamental pulsators, and applying clipping at 
the 2.5ct level (see Table 2, solution 1). They used this 
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Log P 



Fig. 6. Final PL-relation in W for FU mode pulsators in 
the LMC when the depth effect has been taken out. The 
drawn line is the best fit. 




0.5 1 1.5 

Log P 

Fig. 7. As Fig. 6 for FO mode pulsators in the LMC. 



cut-off in period as there is a change of slope observed at 
shorter periods (Bauer et al. 1999). Using the same criteria 
I derive almost identical results (solution 2). A solution 
with clipping at the 4<r level is also included (solution 
3). As mentioned by U99a, and explicitly confirmed here 
(solution 4), the slope for the sample with logP < 0.4 is 
significantly steeper, confirming the result of Bauer et al. 
(1999). 

Table 2 gives the results for FO pulsators as well, in- 
cluding no cut-off in period (entry 5 as obtained by U99a, 
and solution 6 obtained in this paper), and using a cut-off 
at logP = 0.25, for reasons explained before (solutions 7 
and 8). The slope for the shorter periods is significantly 
steeper, as for the FU pulsators. The slopes for the FO 
pulsators are nearly identical to the ones for the FU pul- 
sators both above and below the cut-off. The dispersion in 
all relations is larger than the corresponding relations for 
the LMC, and this is very likely due to the known depth 
effect of the SMC. 

To discuss the spatial structure of the SMC the FU 
pulsators are considered first in determining 9 and i (so- 
lution 3). An arbitrary distance to the center of the SMC 
of R = 63 kpc is assumed, and the center is assumed to 
be given by the mean of the 464 Cepheids in the sample: 




18 16 14 12 10 8 -0.5 

R.A. (degrees) X ^P c ) 

o " ■ o ■ 




-1 1 -0.5 0.5 

Y' (kpc) X' (kpc) 

Fig. 8. As Fig. 2 for the SMC, FU pulsators and i = 
and 9 = 58°. 

a Q = 13.529, d = -72.951 degrees (J2000). Figure 8 gives 
the relevant projection for 9 = 58 and i = degrees. 

It is clear that the SMC is inclined with the NE part 
closer to us. A fit with a position angle of the line-of- 
nodes of 9 = 238 (± 7), and i = 68 (± 2) degrees is 
shown in Fig. 9. The strong correlation in the bottom left 
panel is an artefact. In the transformation to y' and z' , the 
dominant term is respectively — zosini and +ZQCOsi (see 
the Appendix in Weinberg & Nikolaev, 2000). Because, as 
mentioned before, zq does not correspond to the physical 
distance from the plane, but is dominated by the intrin- 
sic scatter in the PL-relation, y' and z' will be strongly 
correlated. 

The analysis is repeated for the FO pulsators and log P 
>0.25 (solution 7), and the result is a P. A. of the line- 
of- nodes 245 ± 10°, and i — 68 ± 4°, consistent with the 
results for the FU mode pulsators. Figure 10 shows the re- 
sult. The values for the inclination angle and the position 
angle of the line-of-nodes are in good agreement with the 
previous results for Cepheids quoted at the beginnning of 
this section, that were obtained using far fewer Cepheids 
but covering a larger area. 

The PL-relations are fitted taking away the depth ef- 
fect, and the results are listed in Table 2 (solutions 9&10), 
and plotted in Figs. 11-12 for FU and FO pulsators. The 
dispersion is lowered significantly but is still larger than 
that in the PL-relation for the LMC Cepheids, indicating 
a substantial intrinsic depth of the SMC. 

The total depth effect, if the intrinsic scatter in the 
PL-relation were zero, over all stars in the direction of 
the SMC bar is about 20 kpc (Fig. 13), with the depth 
along any given line-of-sight towards the bar being about 
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Table 2. PL-relations in the SMC of the form M = a x logP + b 



solution 


M 


a 


b 


a 


N 


Remarks 


1 


W 


-3.310 ±0.020 


16.387 ± 0.016 


0.125 


463 


U99a; FU; log P > 0.4; 2.5cr clipping 


2 


w 


-3.300 ±0.021 


16.381 ± 0.016 


0.126 


446 


This paper; FU; log P > 0.4; 2.5cr clipping 


3 


w 


-3.286 ± 0.024 


16.366 ± 0.018 


0.145 


464 


This paper; FU; log P > 0.4; 4a clipping 


4 


w 


-3.539 ±0.058 


16.520 ± 0.014 


0.142 


743 


This paper; FU; log P < 0.4; 4a clipping 


5 


w 


-3.567 ±0.023 


15.981 ± 0.006 


0.121 


714 


U99a; FO; logP > -0.2; 2.5cr clipping 


6 


w 


-3.558 ± 0.025 


15.978 ± 0.006 


0.126 


688 


This paper; FO; logP > —0.2; 2.5cr clipping 


7 


w 


-3.284 ±0.091 


15.854 ± 0.037 


0.138 


211 


This paper; FO; logP > 0.25; 3a clipping 


8 


w 


-3.572 ± 0.050 


15.982 ± 0.067 


0.145 


518 


This paper; FO; logP < 0.25; 3a clipping 


9 


w 


-3.328 ± 0.023 


16.395 ± 0.017 


0.132 


464 


FU; log P > 0.4; 4a clipping; planar effect taken out 


10 


w 


-3.360 ± 0.082 


15.881 ± 0.033 


0.123 


211 


FO; logP > 0.25; 3a clipping; planar effect taken out 




16 kpc (see Fig. 8 lower left panel). These values are in 
fact upper limits because of the intrinsic scatter in the PL- 
relation. Previously it was shown that the intrinsic depth 
of the LMC is negligible and that the scatter in the LMC 
PL-relation is intrinsic. If it is assumed that the intrinsic 
scatter in the SMC PL-relation is the same, this would 
imply that the intrinsic depth of the SMC translates into a 
dispersion of about V0.132 2 — 0.072 2 = 0.11 mag meaning 
that the above numbers should be multiplied by 0.7 to 
get the true depth effect. The estimates for the depth of 
the SMC are consistent with those based on independent 
analysis, as quoted in Westerlund (1997, his Sect. 12.2.1). 

5. Infrared PL-relations 

In this section infrared PL-relations in the SMC and LMC 
will be determined, based on the OGLE Cepheids that have 
a counterpart in the 2MASS survey. 




0.5 1 1.5 

Log P 



Fig. 11. Final PL-relation in W for FU mode pulsators 
with logP > 0.4 in the SMC when the depth effect has 
been taken out. The drawn line is the best fit. 



The observed JHK magnitudes are first corrected for 
reddening. The E(B — V) of the respective OGLE field is 
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Fig. 12. As Fig. 11 for FO mode pulsators with logP > 
0.25 in the SMC. 



LO 




Z (kpc) 

Fig. 13. If the intrinsic spread of the Pi-relation were 
zero, this would be the distribution of stars along the 
line-of-sight towards the SMC. Positive- z are towards the 
observer. As the intrinsic spread is about 0.11 mag, the 
horizontal scale should be compressed in such a way that 
the 10 kpc tick corresponds to 7 kpc. 

Table 3. Derived values for the inclination and position 
angle of the line-of-nodes 



Galaxy 


i 


e 


LMC 


18 ±3° 


296 ± 1° 


SMC 


68 ±2° 


238 ± 7° 



used (U99b,c), a selective reddening Ay/E(B — V) of 3.1 
is used, and the extinction curve of Cardelli et al. (1989). 
In particular, At = 0.90P(P - V), Ah = 0.59P(P - V) 
and A K = 0.36P(P - V) are used. Second, the 2MASS 
magnitudes are transformed to the Carter system using 
the average formulae derived in Sect. 3 (that is, -fT(Carter) 
- if (2Mass) = +0.02, PT(Carter) - P(2Mass) = -0.02, 
J(Carter) - J(2Mass) = +0.05). Third, the planar effect 
is taken out using the adopted values in Table 3 based on 
the results obtained from the FU pulsators. 

One issue that is important here is Malmquist bias. 
For short periods only the on average brighter Cepheids 
have been detected by 2MASS because the magnitudes are 
near the detection limit of 2MASS. Not taking this effect 
into account would lead to a derived slope of the PL- 



relation which is too shallow (see Lanoix et al. 1999b). 
The nominal completeness limits are J = 15.8, H = 15.1 
and K = 14.3 mag but at higher galactic latitude the 
2MASS catalog contains accurate detections 0.5-1.0 mag 
fainter than this (Cutri et al. 2000). 

Table 4 contains the derived PL-relations in the dif- 
ferent bands, on the Carter system, and with the planar 
effect taken out. Different cut-offs in logP were tested 
to ensure that the slopes are unbiased. The PL-relations 
are displayed in Figs. 14-25. The dispersions are generally 
slightly larger than quoted in the literature (for exam- 
ple, 0.16, 0.12, 0.12 mag in the PL-relation for FU pul- 
sators in respectively J, H, K [Gieren et al. 1998, using 59 
Cepheids in the LMC]), but this is due to the fact that 
single-epoch infrared photometry is used here, while the 
data that Gieren et al. use is based on intensity-mean mag- 
nitudes. Peak-to-peak amplitudes depend on period and 
wavelength but can be 0.2-0.6 mag in the infrared (Laney 
& Stobie 1986b). The fact that the sample studied here is 
an order of magnitude larger than previously studied still 
allows for very accurate determinations of the slope and 
zero point of the infrared PL-relations. 

Some tests were made to check the influence of the 
adopted reddening. Decreasing E(B — V) in steps of 0.01 
mag will make the respective zero points in JHK fainter 
in steps of respectively 0.009, 0.006 and 0.004 mag. For 
reference, based on the solutions in the P-band, the av- 
erage E(B - V) is 0.148 for the FU pulsators, and 0.156 
for the FO pulsators in the LMC (averaged over all pe- 
riods), and respectively 0.091 and 0.092 for the FU (for 
logP > 0.4) and FO (for logP > 0.3) pulsators in the 
SMC. 

There is some evidence that the selective reddening 
towards the LMC is lower than the adopted value of 3.1 
(Misselt et al. (1999) derive values for R towards LMC 
stars that range between 2.16 ± 0.30 and 3.31 ± 0.20, 
with un unweighted mean of 2.55). Chosing R = 2.6 would 
make the respective zero points in JHK fainter by respec- 
tively 0.037, 0.021 and 0.005 mag in the LMC, and fainter 
by respectively 0.008, 0.008 and 0.006 mag in the SMC. 

6. Discussion 

6.1. The slopes of the PL—relations 

One of the main results of this paper is the derivation 
of the slopes of PL-relation in the Wesenheit-index (W), 
and in JHK, for both FU and FO pulsators. In this sec- 
tion the slopes are intercompared, compared to previous 
determinations, and to predictions of theoretical models. 

The first thing to remark is that in the LMC the slope 
of the PL-relation of the FU and FO pulsators disagree 
at the 1.6-2.9cr level. For all 4 bands the slope of the PL- 
relation of the FO pulsators is steeper than that of the 
FU pulsators. In the SMC, the slope of the PL-relation of 
the FO pulsators is steeper than that of the FU pulsators 
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Table 4. Infrared PL-relations of the form M = a x log P ± b 
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solution 


M 


a 




b 






a 


N 


Remarks 




1 


K 


-3.246 ± 0.036 


16.032 


± 


0.025 





168 


472 


LMC; 


FU; 


log P > 0.4; 3a clipping 


2 


Ho 


-3.236 ± 0.033 


16.048 


± 


0.023 





161 


493 


LMC; 


FU; 


all log P; 3cr clipping 


3 


Jo 


-3.144 ± 0.035 


16.356 


± 


0.025 





173 


490 


LMC; 


FU; 


all logP; 3<r clipping 


4 


K 


-3.381 ± 0.076 


15.533 


± 


0.032 





137 


238 


LMC; 


FO; 


log P > 0.25; 3<t clipping 


5 


Ho 


-3.434 ± 0.060 


15.550 


± 


0.024 





126 


263 


LMC; 


FO; 


all logP; 3a clipping 


6 


Jo 


-3.299 ± 0.066 


15.817 


± 


0.027 





140 


267 


LMC; 


FO; 


all log P; 3a clipping 


7 


Ko 


-3.212 ± 0.033 


16.494 


± 


0.026 





194 


418 


SMC; 


FU; 


log P > 0.4; 3a clipping 


8 


Ho 


-3.160 ±0.032 


16.475 


± 


0.025 





182 


414 


SMC; 


FU; 


logP > 0.4; 3a clipping 


9 


Jo 


-3.037 ±0.034 


16.771 


± 


0.027 





199 


418 


SMC; 


FU; 


log P > 0.4; 3a clipping 


10 


Ko 


-3.102 ± 0.155 


15.937 


± 


0.068 





178 


156 


SMC; 


FO; 


logP > 0.3; 3a clipping 


11 


Ho 


-3.199 ±0.142 


15.990 


± 


0.062 





162 


162 


SMC; 


FO; 


log P > 0.3; 3a clipping 


12 


Jo 


-3.104 ±0.159 


16.292 


± 


0.070 





183 


158 


SMC; 


FO; 


logP > 0.3; 3a clipping 




Log P 



Fig. 14. Final PL-relation in K for FU mode pulsators in 
the LMC when the depth effect has been taken out. The 
drawn line is the best fit. 





Log P 



Fig. 16. As Fig. 14: Final PL-relation in J for FU mode 
pulsators in the LMC. 



Fig. 15. As Fig. 14: Final PL-relation in H for FU mode 
pulsators in the LMC. 



in WJH, but in all 4 bands the difference is less than la, 
largely influenced by the fact that the error in the slope 
of the PL-relation for the SMC FO pulsators is relatively 
large. 

The second thing to remark is that for the FU pul- 
sators the slope for the PL-relation of the LMC and SMC 
Cepheids agree within la in WK. In JH the slope of the 
PL-relation for the LMC pulsators is steeper than that 
for the SMC pulsators at about the 1.1 -2. 2a level. For the 



Fig. 17. As Fig. 14: Final PL-relation in K for FO mode 
pulsators in the LMC. 

FO pulsators the slope of the PL-relation for the LMC 
Cepheids is steeper than that of the SMC Cepheids in all 
4 bands but the difference is 1.5a in HK and less in W J . 

Table 5 collects slopes quoted in the literature both 
from observations and theoretical models, and the fi- 
nal values derived in the present paper. Regarding the 
Wesenheit-index, the only previous determination was by 
Tanvir (1999), who used a different definition for this in- 
dex and so can not be compared directly to the present 
result. The results by U99a are essentially based on the 
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Fig. 18. As Fig. 14: Final PL-relation in H for FO mode 
pulsators in the LMC. 
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Log P 

Fig. 19. As Fig. 14: Final PL-relation in J for FO mode 
pulsators in the LMC. 
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Fig. 20. As Fig. 14: Final PL-relation in K for FU mode 
pulsators in the SMC. 

same sample considered here and hence are very close to 
the present results. 

Regarding the IR relations, the slopes derived in the 
present work are the most accurate available for the mo- 
ment. The slopes agree with in the errors with all previ- 
ous determinations. The previously most accurate slopes 
were by Gieren et al. (1998)0. Also here the agreement is 
excellent. The reliability of the result is enhanced if one 

4 I have calculated the error in the slope using the data they 
kindly provided (their Table 9), and made a division on pe- 
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Log P 

Fig. 21. As Fig. 14: Final PL-relation in H for FU mode 
pulsators in the SMC. 
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Log P 

Fig. 22. As Fig. 14: Final PL-relation in J for FU mode 
pulsators in the SMC. 




0.5 1 1.5 

Log P 

Fig. 23. As Fig. 14: Final PL-relation in K for FO mode 
pulsators in the SMC. 



considers that the Gieren et al. sample contains mainly 
long-period Cepheids while the OGLE sample is dominated 
by short-period Cepheids. 

The agreement with theoretically predicted slopes is 
less satisfactory. For the Wesenheit-index the predicted 
slopes by Bono & Marconi (2000, private communi- 
cataion) are too shallow compared to the observations at 
the 2-4cr level. In J the models by Alibert et al. (1999) 

riod as the sample they consider almost exclusively contains 
Cepheids with long periods, unlike the OGLE sample. 
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Table 6. The relative distance modulus SMC-LMC 



Log P 

Fig. 24. As Fig. 14: Final PL-relation in H for FO mode 





A 




Remarks 


0.51 


± 


0.03 


U99c; W; FU 


0.510 


± 


0.023 


present work; W; FU 


0.522 


± 


0.042 


present work; W; FO 


0.479 


± 


0.044 


present work; K; FU 


0.488 


± 


0.091 


present work; K; FO 


0.465 


± 


0.041 


present work; H; FU 


0.511 


± 


0.081 


present work; H; FO 


0.469 


± 


0.044 


present work; J; FU 


0.534 


± 


0.092 


present work; J; FO 


0.475 


± 


0.022 


mean of the 6 IR determinations 


0.504 


± 


0.015 


mean of all 8 determinations 



pulsators in the SMC. 

Table 7. Metallicity dependence of the absolute mag- 

2 | — ■ — . — ■ — . — , — . — . — ■ — . — , — . — ■ — ■ — ■ — , — . — | nitude between LMC and SMC, AM = M(LMC) 

-M(SMC) at logP = 0.5 for FU pulsators. 



AM, 



AMk Reference 



+0.030 



- -0.093 



+0.020 Bono et al. (1999) 
-0.091 Alibert et al. (1999) 



CD 
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Log P 

Fig. 25. As Fig. 14: Final PL-relation in J for FO mode 
pulsators in the SMC. 

are clearly too steep. In K the agreement with the Bono 
et al. (1999) models is within the errors, but the Alibert 
et al. (1999) models predict slopes that are too steep. 

Bono et al. (1999) are the only ones to predict the 
slopes for FO-pulsators (averaged over Z = 0.02, 0.008 
and 0.004). Their prediction is in good agreement with the 
observed slope in the LMC and within 2a of the observed 
slope in the SMC. 

Clearly, the present results can be used as tight con- 
straint for theoretical modelling as for the first time accu- 
rate slopes in the infrared for both FU and FO pulsators 
are presented. 

6.2. The relative distance SMC-LMC 

These results allow a determination of the relative distance 
SMC to LMC, based on the various PL-relations. Since 
the slopes of the PL-relations in SMC and LMC are not 
exactly the same, the absolute magnitudes are calculated 
for logP = 0.5 when using FU pulsators, and logP = 
0.3 for FO pulsators for each Cloud and then subtracted, 
using the PL-relations where the distance effect has been 
taken out. The results are listed in Table 6. The errors 
are based on the errors in the slope and zero point. It 
seems that the difference in distance modulus is somewhat 



less using the IR PL-relations than from the Wesenheit- 
index. A mean using all available 8 determinations yields 
Asmc-lmc = 0.50+0.02 mag. This is in good agreement 
with the original determination in U99c based primarily 
on the Wesenheit index for the FU pulsators, and the 0.50 
± 0.03 determined by Cioni et al. (2000b) as the average 
of the difference in J UK and rrtboi of the TRGB (Tip of 
the Red Giant Branch) magnitudes based on DENIS data 
for the two Clouds. 

All this assumes no explicit metallicity correction. As 
recapitulated in GO00 there is no agreement between the- 
ory and empirical evidence, and amongst different theoret- 
ical models, about the extent of a metallicity dependence 
of the PL-relation, if any. In the infrared and Wesenheit- 
index there are no empirical estimates for the metallic- 
ity dependence. Table 7 summarises the predictions of 
AM = M(LMC) -M(SMC) in W, J, K for 2 sets of recent 
models at logP = 0.5. Again the disagreement is obvious 
and at a level of upto 0.1 mag which is not negligible. 

7. The zero point of the Galactic PL-relation 

GO00 discuss in detail the sample of Galactic Cepheids 
with hipparcos parallaxes and the methods and rela- 
tions used to compute the zero point of the Galactic PL- 
relation. Since the slopes in W and K differ from those 
used in GO00, new zero points calculations are necessary. 
They are summarised in Table 8, and have the same for- 
mat as in GO00. It was remarked in GO00, and repeated 
here, that the zero point of both the V and / PL-relations 
for the sample of the 63 stars with infrared data differs by 
0.10 mag from those using the full samples. Until infrared 
date for more Cepheids in HIPPARCOS becomes available 
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Table 5. Slopes of the Pi-relations. 



Slope 


Colour 


System 


Reference 






Emperical results, FU pulsators 


-3.337 ±0.013 


Wesenheit 


LMC 


this work 


-3.411 ±0.036 


Wesenheit 


LMC 


Tanvir (1999); for W = I - 1.45 (V - 7) 


-3.300 ±0.011 


Wesenheit 


LMC 


Udalski et al. (1999a) 


-3.328 ±0.023 


Wesenheit 


SMC 


this work 


-3.310 ±0.020 


Wesenheit 


SMC 


Udalski et al. (1999a) 


-3.144 ±0.035 


J 


LMC 


this work 


-3.129 ±0.052 


J 


LMC 


Gieren et al. (1998) 


-3.147 ±0.065 


J 


LMC 


this work; 54 stars with logP > 0.833 from Gieren et al. (1998) 


-2.99 ±0.59 


J 


LMC 


this work; 5 stars with logP < 0.833 from Gieren et al. (1998) 


-3.31 ±0.10 


J 


LMC 


Madore & Freedman (1991) 


-3.31 ±0.12 


J 


LMC 


Laney & Stobie (1986a; their solution lb) 


-3.037 ± 0.034 


J 


SMC 


this work 


-3.22 ±0.17 


J 


SMC 


Laney & Stobie (1986a; their solution 10b) 


-3.236 ± 0.033 


H 


LMC 


this work 


-3.249 ± 0.044 


H 


LMC 


Gieren et al. (1998) 


-3.246 ± 0.056 


H 


LMC 


this work; 54 stars with logP > 0.833 from Gieren et al. (1998) 


-3.49 ±0.26 


H 


LMC 


this work; 5 stars with logP < 0.833 from Gieren et al. (1998) 


-3.72 ±0.07 


H 


LMC 


Caldwell & Laney (1991) 


-3.40 ±0.10 


H 


LMC 


Laney & Stobie (1986a; their solution 2b) 


-3.37 ±0.10 


H 


LMC 


Madore & Freedman (1991) 


-3.160 ±0.032 


H 


SMC 


this work 


-3.25 ±0.07 


H 


SMC 


Caldwell & Laney (1991) 


-3.36 ±0.15 


H 


SMC 


Laney & Stobie (1986a; their solution lib) 


-3.246 ± 0.036 


K 


LMC 


this work 


-3.267 ± 0.041 


K 


LMC 


Gieren et al. (1998) 


-3.304 ± 0.052 


K 


LMC 


GO00 [54 stars with logP > 0.833 from Gieren et al. (1998)] 


-3.37 ±0.39 


K 


LMC 


GO00 [ 5 stars with logP < 0.833 from Gieren et al. (1998)] 


-3.42 ±0.09 


K 


LMC 


Madore & Freedman (1991) 


-3.41 ±0.10 


K 


LMC 


Laney & Stobie (1986a; their solution 3b) 


-3.212 ± 0.033 


K 


SMC 


this work 


-3.38 ±0.15 


K 


SMC 


Laney & Stobie (1986a; their solution 12b) 






Theoretical results, FU pulsators 


-3.02 ±0.04 


Wesenheit 


0.02 


Bono & Marconi (2000, private communication) 


-3.17 ±0.04 


Wesenheit 


0.008 


Bono & Marconi (2000, private communication) 


-3.21 ±0.05 


Wesenheit 


0.004 


Bono & Marconi (2000, private communication) 


-3.286 


J 


0.01 


Alibert et al. (1999) 


-3.262 


J 


0.004 


Alibert et al. (1999) 


-3.19 ±0.09 


K 


0.008 


Bono et al. (1999) 


-3.395 


K 


0.01 


Alibert et al. (1999) 


-3.27 ±0.09 


K 


0.004 


Bono et al. (1999) 


-3.369 


K 


0.004 


Alibert et al. (1999) 






Emperical and theoretical results, FO pulsators 


-3.381 ± 0.076 


K 


LMC 


this work, FO-pulsators 


-3.102 ± 0.155 


K 


SMC 


this work, FO-pulsators 


-3.44 ±0.05 


K 


all Z 


Bono et al. (1999), FO-pulsators 



to better constrain this zero point, it seems prudent to 
add 0.10 mag to the zero points in K as listed in Table 8. 

As discussed in GO00, when taking into account the 
errors in reddening, period-colour relations, etc., an er- 
ror term of 0.028 for W and 0.011 in K should be added 
in quadrature to the quoted errors, and the zero points 
should also be corrected for Malmquist-bias (±0.01 mag) 
as derived in GO00 from theoretical calculations and nu- 
merical simulations. The final zero points for the 4 val- 



ues listed in Table 8 become respectively —2.711 ± 0.108, 
-2.718 ± 0.108, -2.517 ± 0.169 and -2.545 ± 0.169. 

These zero points can now be combined with the zero 
points of the LMC and SMC PL-relations to yields dis- 
tances of 18.60 ± 0.11 (W) and 18.55 ± 0.17 (K) for the 
LMC, and 19.11 ± 0.11 (W) and 19.04 ± 0.17 (K) for 
the SMC. The distance to the LMC based on the Wesen- 
heit index is identical to the value derived in GO00, which 
was based on the results of Tanvir (1999) for the Wesen- 
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heit index for a much smaller sample of LMC Cepheids. 
To these distances, systematic uncertainties must be con- 
sidered due to the fact that the slope of the galactic PL- 
relation may not be the same as in the MCs and the uncer- 
tain metallicity effect. In GO00 these unceratinties were 
estimated to be (±0.08 slope) (^0.15 metallicity) in W, and 
(±0.03 slope) (±0.06 metallicity) in K. 

I do not want to enter here into a discussion about the 
"long" and "short" distance scale to the LMC (see Gibson 
1999 for an overview). The latest "short" distance modu- 
lus to the LMC based on the Red Clump method stands 
at 18.27 ± 0.07 (Popowski 2000) or 18.24 ± 0.08 (Udal- 
ski 2000). That to the SMC at 18.77 ± 0.08 (Popowski 
2000). The distance to the LMC and SMC derived here 
are "long" and in excellent agreement with those derived 
recently by Cioni et al. (2000b) using the TRGB-method 
from denis data of respectively 18.55 ± 0.04 (formal) 
± 0.08 (systematic) and 18.99 ± 0.03 (formal) ± 0.08 
(systematic), and with the distance to the SMC of 19.05 
± 0.017 (standard error) ± 0.043 (systematic error) us- 
ing FU/FO and FO/SO (Second Overtone) double-mode 
Cepheids (Kovacs 2000) and 19.11 ± 0.08 using pure SO 
and FO/SO double-mode Cepheids (Bono et al. 2000). 
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Appendix A 

In this Appendix a comparison is presented between the 
DENIS / and OGLE / and denis JK and 2MASS JK for 
the respectively 173 and 141 LMC and SMC Cepheids in 
common. This has its limitations since the comparison is 
done using variable stars. Yet it may still be of interest 
to other workers, as such a comparison has not been done 
yet. The variability would increase the spread in any cor- 
relation but should not wash out any colour terms, if they 
exist. 

The observed UK magnitudes are first corrected for 
reddening. The E(B — V) of the respective OGLE field is 
used (U99b,c), a selective reddening A V /E(B -V) of 3.1 
is used, and the extinction curve of Cardelli et al. (1989). 



In particular, A l = 1.87 E(B - V), Aj = 0.90P(P - V) 
and Ak = 0.36E(B—V) are used. Since the effective wave- 
lengths of the denis and OGLE I and the denis and 2MASS 
JK are very similar, no distinction between the respective 
filterscts is made. At the present level of accuracy this ef- 
fect of differential reddening between the slightly different 
filters should be entirely negligible. 

Figure Al shows the differences between 2 MASS and 
DENIS J and K plotted versus 2MASS (J - K) Q , and 
OGLE and denis I a and plotted versus DENIS (/ - J) - In 
the left panel, all stars are plotted, in the right panel a 
sub-sample is plotted using a criterium on the total er- 
ror. Linear least-square fits were made and the results are 
listed in Table Al. 

From this data, there is no evidence for a colour-term 
in I and J and the difference between OGLE and denis 
Io, and 2MASS and denis Jo is less than 0.10 magnitude. 
In K there seems to be a color term at the 3a level when 
all stars are used, which still is present at the 2a level 
when a subsample of stars with the smaller photometric 
errors is being used. The average difference if (2Mass) 
— Ao(Denis) is +0.20 magnitudes, which is comparable 
to the observed scatter, which is an upper limit to the 
intrinsic scatter as two single-epoch data are being com- 
pared. 

The conclusion is drawn from the present dataset that 
there is no evidence for colour terms when comparing 
OGLE and denis Jo, and 2MASS and DENIS J , and that 
any difference between the systems is less than 0.10 mag- 
nitude. Regarding 2MASS and denis K 0} there is evi- 
dence for a colour term, which amounts to a difference 
AT (2Mass) — AT (Denis) = +0.20 magnitudes for a typi- 
cal colour of (J — K) = 0.4. 
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Table 8. Values for the zero point for Galactic PL- relations 



Colour 


N 
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Assumed 


Total 


Remarks 








Slope 


Weight 
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W 


191 


-2.728± 0.104 


-3.328 


5331.0 


All stars 


K 
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K 
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Table Al. Comparing photometric systems. Fits of the form y — a x x ± b 
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a 




X 


b 


a 


N 


K (2Mass)-K (Denis) 
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(J- 
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